Spin-orbital coupling in a triplet superconductor-ferromagnet junction 
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We study a novel type of coupling between spin and orbital degrees of freedom which appears 
at triplet superconductor-ferromagnet interfaces. Using a self-consistent spatially-dependent mean- 
field theory, we show that increasing the angle between the ferromagnetic moment and the triplet 
vector order parameter enhances or suppresses the p-wave gap close to the interface, according as the 
gap antinodes are parallel or perpendicular to the boundary, respectively. The associated change in 
condensation energy establishes an orbitally-dependent preferred orientation for the magnetization. 
When both gap components are present, as in a chiral superconductor, we observe a first-order 
transition between different moment orientations as a function of the exchange field strength. 

PACS numbers: 74.45. +c, 74.20. Rp, 74.50.+r 



Introduction. The singlet superconductor (SSC) and 
ferromagnet (FM) phases are fundamentally incompati- 
ble, as the exchange field of the FM acts to align the anti- 
parallel spins of the electrons in singlet Cooper pairs, thus 
destroying the superconductivity Due to this pair- 
breaking effect, homogeneous coexistence of SSC and FM 
is very rare. On the other hand, SSC-FM interfaces can 
be readily fabricated in artificial layered heterostructures, 
and the study of these devices has attracted intense at- 
tention 043 • The pair- breaking effect is central to the 
understanding of these systems, e.g. it causes the spatial 
oscillation of the SSC correlations in the barrier of a ferro- 
magnetic Josephson junction, which is responsible for the 
famed 0-7T transition 0, 01] . It is well known that the FM 
suppresses the SSC gap close to the interface 0,0,0], and 
can induce a magnetization in the SSC jij]. Conversely, 
the SSC may influence the FM in order to minimize pair 
breaking, by suppressing the uniform magnetization near 
to the interface [4| , or causing the spontaneous formation 
of domains in a thin FM layer so as to give zero averaged 
magnetization on the scale of the coherence length [5[ ■ 



The coexistence of FM and triplet superconductor 
(TSC) states is more favorable, as the exchange field is 
only pair breaking when it is perpendicular to the Cooper 
pair spins. The physics of TSC-FM devices is therefore 
richer than their singlet counterparts, as the orientation 
of the FM moment relative to the TSC vector order pa- 
rameter is now a crucial variable. This is predicted to 
control the nature of the proximity effect in TSC-FM bi- 
layers H and the sign of the current in TSC-FM-TSC 
Josephson junctions [9(. In addition to this orientation- 
dependent pair breaking, the interface with the FM also 
induces distinctly new physics within the TSC [Io|, EI 



spin-flip reflection processes scatter the triplet Cooper 
pairs between the spin f and 4- condensates of the TSC, 
setting up a Josephson-like coupling between them. The 
resulting "spin Josephson effect" is manifested as a spon- 



taneous s pin current in the TSC normal to the TSC-FM 
interface [lfj, and is expected to contribute to the inter- 
face free energy The study of TSC-FM interfaces 
is timely, as the recent preparation [12| of superconduct- 
ing thin films of Sr 2 Ru04 [l3| opens the way to TSC 
heterostructures. Furthermore, the proposed appearance 
of Majorana fermions at TSC-FM interfaces in quantum 
wires motivates a deeper understanding of the interplay 
between FM and TSC " 



14] 



In this letter we demonstrate that, contrary to what is 
expected from the pair- breaking effect, it is the orbital 
and not the spin component of the TSC order parame- 
ter which dictates the energetically-favorable orientation 
of the FM moment in a TSC-FM heterostructurc. This 
reveals the emergence of a distinct form of spin- orbital 
coupling at TSC-FM interfaces, which originates from the 
spin-flip scattering of triplet Cooper pairs. We study a 
lattice model of a TSC-FM hybrid using a self-consistent 
Bogoliubov-de Gennes theory 0, [HI , which accounts for 
the pair-breaking and spin-flip scattering processes on 
an equal footing. For a unitary p-wave TSC, we show 
that the latter stabilizes the FM moment either paral- 
lel or perpendicular to the spin orientation of the Cooper 
pairs, depending on the alignment of the p-wave gap with 
respect to the interface. For the chiral p x + ip y state, in- 
creasing the strength of the FM exchange leads to a first 
order transition between these two configurations. 

The model. We study a finite-size square lattice 
model of the planar TSC-FM heterostructure shown 
in Fig. [lja). The lattice size is L x L, with periodic 
boundary conditions imposed along the direction par- 
allel to the interface. Indicating each site by a vector 
i = (i x ,i v ), with i x and i y being integers ranging from 
—L/2 to L/2, we write the Hamiltonian 
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FIG. 1. (a) Schematic diagram of the two-dimensional TSC- 
FM junction. The FM region is located at x < 0, while the 
TSC is realized for x > 0. The magnetization M of the FM 
(big red arrow) forms an angle <j> with the d-vector of the TSC 
(thin blue arrow), which defines the z axis. We study TSC 
states with p x , p y and p x + ip y symmetry, (b) Evolution of 
the bulk FM magnetization and the majority and minority 
spin concentrations as a function of the exchange field h. 
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where c\ a is the annihilation operator of an electron with 
spin a at the site i, (i, j) indicates nearest-neighbor sites, 
fi is the chemical potential, and Si = J2 S s' c t s a s,s' c is' 
is the spin density at site i. The lattice is divided into 
three regions: the FM subsystem for i x < 0, the TSC 
subsystem for i x > 0, and the interface at i x = 0. For 
simplicity, we assume that the hopping matrix element t 
is the same across the lattice; relaxing this assumption is 
not expected to qualitatively alter our results. All energy 
scales are expressed in units of t. 

A nearest-neighbor attractive interaction —V (V > 0) 
is present only on the TSC side of the junction. By tun- 
ing the electron density and pairing strength, at mean- 
field level this gives a TSC state with d-vector paral- 
lel to the z-axis, see Fig. QJa). We consider three or- 
bital symmetries: p x , p y , and p x + ip y . The interaction 
term in H is decoupled within the Hartree-Fock approx- 
imation by introducing the pairing amplitude on a bond 
Ay = (cifCji). The pairing amplitudes are used to con- 
struct the superconducting order parameter with the de- 
sired symmetry in the spin-triplet channel, which is then 
determined self-consistently jg, [l5| . The p x + ip y state 
(relevant for S^RuCi) is energetically most favorable as 
it fully gaps the 2D Fermi surface; anisotropic pairing in- 
teractions with V x = —V, V y = and vice versa are used 
to model the p x and p y states, respectively. We choose 
the bulk TSC critical temperature to be ksT c = 0.18. 
The FM subsystem is modeled by the constant exchange 
field h, which forms the angle (f> with respect to the 
direction of the d-vector as shown in Fig. [Ha). Since 
the TSC state is invariant under spin rotations about 
the d-vector, h is restricted to lie in the x-z plane, i.e. 
h = h(sm(4>), 0, cos{<p)). The relation between the mag- 
netization M and h is shown in Fig. (Ub). 

From the resulting mean- field Hamiltonian Hmf one 
can determine the condensation energy of the TSC 




FIG. 2. Zero temperature pairing amplitude scaled to its bulk 
value as a function of the distance i x from the interface at h = 
1.5 and at different angles tf> of the d-M misalignment. The 
spin-triplet orbital symmetry is of p x (a), p y (b), and chiral 
type (i.e. with a real p x (c) and an imaginary p y component 
(d) ), respectively. 



as well as the total Gibbs energy F of the junction: 
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where (3 = (/c^T)" 1 . Both F and Ea depend on the ex- 
change field and the angle (f> through the self-consistent 
pairing amplitude and magnetization. The results pre- 
sented here were obtained using L = 120; a larger lattice 
does not qualitatively change our conclusions. 

Pairing amplitude. In Fig. [5]we present representative 
examples of the pairing amplitude profile near the inter- 
face for h = 1.5 and several different values of < < \. 
Distinct trends are evident both in the FM and TSC sides 
of the junction as the exchange field is rotated from a par- 
allel {(f> = 0) to a perpendicular (</>=?•) orientation with 
respect to the d-vector. Independent of the orbital sym- 
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metry, the proximity effect in the FM smoothly evolves 
from a monotonous decay to a damped oscillating be- 
haviour as 4> is decreased from ^ to 0. The oscillating 
behaviour is similar to that observed in a SSC-FM junc- 
tion and here it is also due to pair breaking, specifically 
the coupling of the z-component of the exchange field to 
the in-plane spin of the triplet Cooper pairs. 

In contrast, the pairing amplitude on the TSC side 
of the interface strongly depends upon both the angle 
</> and the orbital symmetry of the TSC. Specifically, as 
can be seen in Fig. [Ha), a TSC with p x orbital symmetry 
shows a reduction of the pairing amplitude close to the 
interface as the exchange field is tilted from parallel to 
perpendicular with respect to the d- vector. For &p y TSC 
we observe the opposite behaviour, although the effect is 
less pronounced [see Fig. |5{b)]. The chiral p x + ip y TSC 
evidences both trends: decreasing <f> from ^ to enhances 
the real (p x ) part of the gap [Fig.[2fc)], but suppresses the 
imaginary (p y ) part [Fig. HJd)]. Interestingly, for a given 
exchange field magnitude h, the variation of each gap 
component of the p x + ip y state with <p is much stronger 
than for the purely p x or p y state. 

The pair-breaking effect on both the FM and the TSC 
side, which is the same for all orbital symmetries, cannot 
explain the different (^-dependence of the p y and p x gap 
profiles. This instead can be understood in terms of the 
spin-flip reflection of triplet Cooper pairs at the interface 
with the FM, which is also crucial for the spin Josephson 
effect 

0,1 

1 1| . In such a scattering process, an incident 
Cooper pair with spin a mutually perpendicular to d and 
M acquires the spin- and orbital-dependent phase shift 
ir — 2a<p + AO: the first two terms are due to the spin- flip, 
while the last is due to the phase change of the TSC gap 
upon specular reflection. Here AO = it (0) for the p x (p y ) 
state, while AO depends on the angle of incidence for the 
Px + ipy gap. It is well known that the gap is suppressed 
at interfaces where reflected Cooper pairs undergo a non- 
trivial phase shift [l^]; by choosing <j) so that the spin- 
flip reflected Cooper pairs have a 27m phase shift, we 
hence expect to maximize the order parameter near the 
interface. Due to the different orbital phase shifts AO, 
this occurs at 4> = (■§) for thep x (p y ) pairing amplitude, 
in agreement with the numerics. This interplay of spin 
and orbital degrees of freedom reveals a novel type of 
spin-orbital coupling at the TSC-FM interface. 

Stable moment orientation. The condensation energy 
Ea in the TSC depends upon the angle 4> through the 
variation of the pairing amplitude near to the interface. 
In Fig. EJa) we plot Ea as a function of <f> for several 
typical cases. For the p x and p y gaps the energy is in- 
deed maximized for the exchange field orientation which 
maximizes the gap amplitude. The p x + ip y case is more 
complicated, since here the p x and p y gap components 
show opposite dependence upon (f>. We find that the 
maximum in Ea shifts from cf) = ^ to (f> = with in- 
creasing exchange field strength. That is, for a weak FM 
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FIG. 3. (Color online) (a) Dependence of the condensation 
energy Ea on the angle 4>. (b-d) Dependence of the Gibbs 
energy F on (f> for various fixed h and for (a) p x , (b) p y and (c) 
Px + ipy orbital symmetries of the spin-triplet superconductor. 
-EA.min and Fmin are the minimum amplitudes of the related 
energies. All panels are for the temperature T — 0.05. 



the p y component dominates the physics, while at strong 
polarizations the p x gap is most important. 

The ^-dependent gain in condensation energy implies 
a favored orientation of the FM moment. This is de- 
termined from the Gibbs free energy F of the junction, 
which we plot in Fig. [3Jb)-(d) for the three orbital sym- 
metries as a function of </>. For the p x orbital sym- 
metry, the profile exhibits a single minimum at </> = 
and a maximum at 4> = 7r/2, and vice versa for the p y 
TSC. In the upper panel of Fig. [4] we plot the difference 
F(cp = 0) — F((f) = ^) between these points, which shows 
that the depth of the minimum monotonically increases 
with increasing exchange field h. The stable magnetic 
orientation is therefore parallel (perpendicular) to the d- 
vector if the antinodes of the p-wave TSC gap are per- 
pendicular (parallel) to the interface, see the sketch in 
the bottom panel of Fig. |U 

The Gibbs free energy for the p x + ip y symmetry shows 
a competition between the tendency of the real and imag- 
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FIG. 4. (Color online) Top panel: Gibbs free energy differ- 
ence between the parallel ((/> = 0) and perpendicular (</> = -|) 
configurations of the moment as a function of the exchange 
field strength h at temperature T = 0.05 for the p x , p y , and 
Px + ip y TSCs. Bottom panel: sketch of the most favorable 
magnetic (red arrow) configurations with respect to the or- 
bital symmetry and the d-vector (blue arrow) of the TSC. 



inary gap components to favour an FM with magnetiza- 
tion parallel and perpendicular to the d-vector, respec- 
tively. This interplay produces minima in F at both 
<j> = and 0=7^, with the exchange field strength select- 
ing the stable state. As expected from the condensation 
energy, the <fi — j state is the global minimum up to a 
critical exchange field h cr ^ (« 1.74 at T = 0.05); the the 
= state is stable up to a higher field h cr 2 in the ex- 
treme half-metal regime (« 2.8 at T = 0.05, not shown 
in Fig. h crt i and h cr ^, merge together as the temper- 
ature is increased, so that only the </> = § configuration 
is stable sufficiently close to T c . We note that the com- 
petition between the two gap components is expected to 
be dependent upon the details of the system. 

Our results imply a close connection between the spin- 
orbital coupling and the spin Joscphson effect in the 
TSC-FM junction [13, EH- Specifically, the spin current 
in the Josephson effect is generated by a ^-dependent 
Helmholtz free energy; for a non-self-consistent model of 
a thin FM layer on a bulk TSC, the free energy due to 
the interface electronic reconstruction has similar form 
to that obtained here, and gives the same stable magne- 



tization configurations for the p x and p y states [11 



Experimental considerations. Although the energy dif- 
ference between the (f> = and ^ states shown in Fig. 0] 
is very small, this results from averaging what is essen- 
tially an interface effect over the entire lattice; the energy 
gain per interface unit cell is L times larger, and gives an 



anisotropy energy on the order of ~ 0.01/cbT c for the mi- 
croscopic parameters chosen here. In the case of a thin 
FM layer of width a few unit cells, this is comparable 
to free energy change due to electronic reconstruction 



within the FM itself [11|. The magnetic anisotropy in- 
duced by the coupling to the TSC could be observed by 
ferromagnetic resonance (FMR) measurements: for an 
exchange field h = 0.5 in the FM, and choosing Sr2Ru04 
(T c = 1.5-fiT) for the bulk TSC, we estimate a precession 
frequency of ~ 5 cos(a) GHz. Since there is no spin- 
orbital coupling at SSC-FM interfaces, the observation 
of this precession would strongly indicate a TSC state in 
the superconductor. For a thicker FM layer, we speculate 
that the spin-orbital coupling could modify the magne- 
tization profile near the interface, effectively creating a 
spin-active boundary layer, which can greatly influence 
the proximity effect in the FM Q . 

For the purposes of simplicity we have neglected com- 
plicating factors such as intrinsic anisotropy of the FM, 
demagnetization effects, spin-orbit coupling, etc. Includ- 
ing these requires a more elaborate description, but does 
not change the basic physics of the TSC-FM interface 
discussed here. Some of these effects may even be turned 
to our advantage, e.g. to observe the appearance of the 
anisotropy axis due to the spin-orbital coupling, it would 
be favorable to select an FM with biaxial easy axes ori- 
ented perpendicular and parallel to d. 

Conclusions. In this letter we have evidenced a novel 
type of spin-orbital coupling which appears at TSC-FM 
interfaces. Specifically, an FM moment parallel (perpen- 
dicular) to the TSC d- vector maximizes the magnitude of 
a p-wave order parameter with antinodes perpendicular 
(parallel) to the interface. The origin of this spin-orbital 
coupling is ascribed to the spin-flip reflection of Cooper 
pairs at the TSC-FM interface. We have determined 
the lowest-energy magnetic configuration by calculating 
the Gibbs free energy. Our method allows us to treat 
many different interface physics on an equal footing, e.g. 
the proximity effect, the interfacial electronic reconstruc- 
tion and the superconducting order parameter variation. 
We show that the latter is most important, due to the 
change in the condensation energy of the TSC. For the 
chiral TSC, where both types of p- wave gaps are present, 
the stable moment orientation switches between the per- 
pendicular and parallel configurations with increasing ex- 
change field. The coupling causes the appearance of an 
anisotropy axis in the FM which could be observed in 
FMR measurements, and can act as a test of the orbital 
and spin pairing state of the TSC. 
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